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^ ! 1 Introduction 

^ In the last 20 years we have accumulated an enormous amount of data on 

elementary particles and their interactions. This data serves two purposes: to fix 
the phenomenological parameters of the Standard Model [SM] and to verify that the 
SM is an excellent description of nature. It is our goal to understand the origin of 
these many arbitrary parameters. In this talk we consider a supersymmetric [SUSY] 



a: 



D | SO(10) grand unified theory [GUT]. We present a straightforward procedure, 

incorporating a general operator analysis, which allows us to use low energy data to 
determine the fermionic sector of the theory at the GUT scale 1 . In what follows we 
first review the status quo, i.e. the low energy data [LED]. We then review the 
evidence for SUSY GUTs, and focus on the virtues of the particular group SO(10). 
Following a discussion of our dynamical principles, we present a general operator 
analysis for C } -ermion ,i- e - the fermionic sector of the GUT theory. We argue that all 
fermion masses and mixing angles can be described with a minimum of 5 arbitrary 
parameters in the Yukawa sector of the theory at the GUT scale, Mq- Including the 
parameter tan/3, the ratio of Higgs vevs present in any SUSY theory, we thus obtain 
a 6 parameter description of fermion masses and mixing angles, leading to 8 
predictions. Finally, we present preliminary results. These preliminary results are 
encouraging and eminently testable. Notwithstanding, the power of our analysis is 
in the paradigm whereby the use of LED, hand in hand with symmetry conditions, 
allows us determine the fermion sector of the theory at Mq. 
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2 Status Quo 



There are 18 arbitrary parameters in the SM with 13 of these in the fermion sector of 
the theory.Q Let us briefly review the state of our knowledge of these 18 parameters 
since they play a central role in what follows. In Table 1 we list the 18 parameters of 
the SM along with their experimental and/or theoretical uncertainties. 



Table 1. The 18 parameters of the Standard Model. 



Parameter 


Uncertainty 


Comments 


a, sin 2 6 W 

a s (M z ) 
m e ,m^,m T 
m c ,m b 
m u /m d ,m s /m d 
(m u + m d )/2 

\Vcd\ « \V US \ 

\V c b\, \V u b/V c b\ 


< 1/2 % 
~ 10% 

< 1/3% 

~ 4% 
? 

7 

~ 1.5% 
~ 20% 


high accuracy 
less certain - scale dependent 

high accuracy 
less certain 
chiral Lagrangian - ambiguous 

QCD sum rules 

fair accuracy 

poorly known 


m t ,m H 

m t ,m H 
m t 


m - 150 + 19 + 15 
m t - 150 _ 24 _ 2Q 

m - 131 + 23 + 5 
m - 131 _ 2g _ 5 

m t > 108 GeV 


60 < m H < 1000 in SM 

60 < m H < 150 in MSSM 
Fermilab data 


M z 

J (Jarlskog invariant) 


< .025% 
~ 30% 


high accuracy 
uncertainty in 



Some of these parameters are known to high accuracy, these include a, sin 2 9w, m e , 
m^jrrir, Mz- A few are known with fair accuracy, \V cd \,m c ,rrib. Finally, the follow- 
ing are poorly known: a s (M z ),m u ,md,m s , \V c b\, \V u b/V c b\ ,m t ,rriH and the Jarlskog 
invariant measure of CP violation, J. 

Our knowledge of a s {Mz) suffers predominantly from theoretical uncertainties 
associated with renormalization scale ambiguities. The experimentally allowed range 
is now a s (Mz) = .118±.007 3 , where the errors probably underestimate the theoretical 
uncertainties. 

Light quark masses are determined using chiral Lagrangians and QCD sum rules. 
In the chiral Lagrangian approach, light fermion masses are characterized by their 
transformation properties under the SU(3)l x SU(3)r chiral symmetry. The mass 
matrix, M, transforms as (3,3). However, as shown by Kaplan and Manohar 4 , 

*That is if we assume that neutrinos are massless. Incorporating neutrino masses leads to an 
additional 9 parameters, 3 masses and 6 mixing angles. 



M' = M + Kdet(M)M^~ l with k an arbitrary constant, has the same transformation 
property. This leads to a large theoretical uncertainty in light quark masses. 

The weak mixing angles \V u b\ and \V c b\ suffer from both large theoretical and/or 
experimental uncertainties. The theoretical uncertainties are evidenced by the dif- 
ferent model dependent calculations of B decay. It is hoped that the application of 
heavy quark effective field theory to the exclusive semileptonic B decays will reduce 
these uncertainties, but this requires better statistics. The latest experimental result 
from B — > D*lv is 5 \V c b\ = .050 ± .008 ± .007 where the first error is from statistics 
and the second is from extrapolation uncertainties. The ratio IV^/V^I, determined 
from inclusive B decay, has large model dependences. Recently the experimental re- 
sults from CLEO II changed significantly from previous measurements by both Argus 
and CLEO. The latest results give 5 .038 < |V^&/Vc6| < .097. This is a factor of two 
smaller than previous results and the change was significantly larger than any of the 
experimental errors for any particular model. We will just have to wait and see if this 
settles down with more data. 

The top and Higgs masses are only known from the radiative effects they have on 
electroweak parameters. The latest data, analyzed by Langacker 2 , gives the results 
presented in Table 1, where MSSM denotes the minimal supersymmetric standard 
model. The difference between the SM and MSSM result is the requirement of a 
lighter Higgs boson in the MSSM. 

Finally, J suffers from strong interaction uncertainties associated with the so- 
called bag constant, Bk- The value of JxBk can be derived given the experimental 
value for €k, all quark masses, and the magnitudes of all Kobayashi-Maskawa [KM] 
elements. Thus, the uncertainty in J is limited by the uncertainty in Bk which is of 
order 30%. 

I have briefly reviewed the status quo, since some of these parameters will be 
used as imputs to fix the fundamental parameters in the fermion mass matrices and 
the others test the theory. For example, the biggest uncertainties in our prediction 
of the top mass comes from the uncertainties in and a s (Mz)- Reducing these 
uncertainties would make this prediction much tighter. 

3 SUSY GUTs and the virtues of SO(10) 

• SUSY GUT 

Given two parameters the fine structure constant, ag, at the GUT scale, Mq, 
one determines the three low energy parameters, c^sin 2 ^ and a s , all evaluated at 
some renormalization scale,//, which for convenience we might choose to be Mz- In 
actuality, we use the two best determined parameters, a and sin 2 9w, to fix a G ~ 1/25, 
M G ~ !0 16 GeV and predict (in a SUSY GUT) a s (M z ) = .125 ± .002 ± .009 for a 
central value of m t = 138GeV. 2 []The errors take into account uncertainties in m t 

§In a non-SUSY GUT the prediction is a s (Mz) ~ -07 2 , which is inconsistent with the data. 



and estimates of threshold corrections at both M G and the weak scale. This result is 
in remarkably good agreement with the data. For a heavier top, the central value for 
a s decreases by several percent. This prediction assumes a supersymmetric desert, 
i.e. the only threshold between M z and M G occurs below ~1 TeV due to the new 
spectrum of states encountered in the MSSM. 

In any GUT, the number of fundamental parameters in the gauge sector of the 
theory decreases by one and only in a SUSY GUT does the resulting prediction agree 
with the low energy data [LED]. The powerful assumption of a SUSY desert has an 
equally important consequence. The LED becomes a window into physics at the GUT 
scale, i.e. measurements at the weak scale gives us information about the physics at 
M G . 

• Virtues of SO (10) 

A single family of fermions fits into one irreducible representation — i.e. 16« D 
{ui,di,ei,Ui} with i = 1,2,3 labelling the three families. We take the 3rd 
family to be the top family. 

The two Higgs doublets required in the MSSM fit into one irreducible representation 
- i.e. 10 D {H, H', H 3 , H^} where H, H' are weak doublets necessary for weak 
symmetry breaking and giving masses to quarks and leptons and H 3 , are 
color triplet Higgs which must get mass of order M G to avoid rapid proton 
decay. We will return to this point shortly. 

In order to accomplish the GUT scale symmetry breaking we must have additional 
representations including {45, 16, 16, • • •}, where, for example, a 16 and 16 vev can 
break SO(10) to SU(5) and the 45 vev can then break SU(5) to SU(3) x SU(2) x 
U(l). This may happen at the same scale, M G << M P , or at two separate scales, 
where the first occurs at a scale t>i such that v 5 = M G « t> 10 << M P . 

We have not considered other possible representations which may be relevant for 
GUT symmetry breaking, such as 54, 126, etc. We shall now assume that only the 
45 plays a crucial role in the generation of fermion masses. It is thus necessary to 
elaborate the possible directions the 45 vev may point in the two dimensional space 
of U(l) subgroups of SO(10) which commute with SU(3) x SU(2) x U(l). Although 
there are only two orthogonal directions in this space, we nevertheless consider the 
following 4 possible symmetry breaking vevs — 

(45) x = Ul0 e ia *X 

(45) y = v 5 e iaY Y (1) 

(45) B _ L = v 5 e ta (*-v(B-L) 

(45} T3R = v 5 e ia ^T 3R (2) 

where we have explicitly represented them as two groups of two orthogonal vevs. We 
consider all four since two of these vevs (one in each group) are well motivated. 



X, in Eq. (1), is the U(l) direction which leaves SU(5) invariant. This is why we 
have taken the magnitude of the vev to be v\o, whereas all the others are taken to be 
v 5 since they do not commute with SU(5). 

B — L, in Eq. (2), just measures baryon number minus lepton number. It can 
play a crucial role in splitting the weak doublet and color triplet Higgs multiplets, i.e. 
solving the hierarchy problem. The Higgs doublets carry zero B-L whereas the triplets 
have non-zero B-L. Thus if the Higgs in the 10 gets mass by coupling to this 45, only 
the color triplets will acquire mass at the scale M G . Hence, this vev is expected to be 
a necessary ingredient in any complete SO(10) model which also solves the hierarchy 
problem. 

4 Dynamic Principles 

Let us now discuss the dynamical principles which guide us towards a theory of 
fermion masses. 

0. At zeroth order, we work in the context of a SUSY GUT with the MSSM below 

M G . 

1. We use SO(10) as the GUT symmetry with three families of fermions {16j % = 

1,2,3} and the minimal electroweak Higgs content in one 10. Using SO(10) 
symmetry relations allows us to reduce the number of fundamental parameters. 

2. We will assume that there are family symmetries which enforce zeros of the mass 

matrix, although we will not specify these symmetries at this time. As we 
will make clear shortly, these symmetries will be realized at the level of the 
fundamental theory defined at M P . 

3. Only the third generation obtains mass via a dimension 4 operator. The fermionic 

sector of the Lagrangian thus contains the term Cf D A 33 = A 16 3 10 16 3 . 
This term gives mass to t, b and r. It results in the symmetry relation — A t = 
\ b = \ T = A at Mq. This relation has been studied before by Ananthanarayan, 
Lazarides and Shafi [ALS] 6 and using and m T as input it leads to reasonable 
results for m t and tan j3. 

4. All other masses come from operators with dimension > 4. As a consequence, 

the family hierarchy will be related to the ratio of scales above M G . We will 
show shortly how to understand the higher dimension operators in terms of 
an effective field theory at M G , obtained by integrating out states with mass 
> M G . 

5. [Predictivity requirement] We demand the minimal set of effective fermion 

mass operators at M G consistent with the LED . 



Figure 1: Tree diagram which leads to effective operator when massive states are 
integrated out. 



Let us now consider the general operator basis for fermion masses. Let 
C-f ermion include operators of the form 

Oij = 16i (• • •)„ 10 (• • -) m 16j (3) 

where 

Mq 45 fc+ i---45 ra 
[ ' ' ' h M l P 45 n x ~ l 

and the 45 vevs in the numerator can be in any of the 4 directions, X, Y, B — L, T 3R 
discussed earlier. 

We said that such operators can be considered as the result of integrating out 
states with mass > Mq. For example, you can convince yourself that an operator 
of the form O22 = I62 10 45g 7i-; M ° I62 is generated by the tree graph of Figure 1, 
assuming v±q » V5 ~ Mq. Note it is at this level in the fundamental theory at Mp 
that additional family symmetries are needed to enforce zeros in the mass matrix. It is 
also trivial to evaluate the Clebsch- Gordon coefficients associated with any particular 
operator since the matrices X, Y, B — L, T 3R are diagonal. Their eigenvalues on the 
fermion states are given in Table 2. 



5 Operator analysis 

Our goal is to find the minimal set of fermion mass operators consistent with the 
LED. With any given operator set one can evaluate the fermion mass matrices for up 
and down quarks and charged leptons. One obtains relations between mixing angles 
and ratios of fermion masses which can be compared with the data. It is easy to show, 
however, without any detailed calculations that the minimal operator set consistent 
with the LED is given by 

£ fermion D 33 + 23 + 22 + O12 "22" texture 

or (4) 
O33 + 23 + 23 + O12 "23"' texture 

It is clear that at least 3 operators are needed to give non-vanishing mass to all 
charged fermions, i.e. det(m a ) 7^ for a = u, d, e. That the operators must be in 



Table 2. Quantum numbers of the 
four 45 vevs on fermion states. 
Note, if u denotes a left-handed 
up quark, then u denotes 
a left-handed charge conjugate 
up quark. 





X 


Y 


B-L 


T3R 


u 


1 


1/3 


1 





u 


1 


-4/3 


-1 


-1/2 


d 


1 


1/3 


1 





d 


-3 


2/3 


-1 


1/2 


e 


-3 


-1 


-3 





e 


1 


2 


3 


1/2 


V 


-3 


-1 


-3 





V 


5 





3 


-1/2 



the [33, 23 and 12] slots is not as obvious but is not difficult to show. It is then easy 
to show that 4 operators are required in order to have CP violation. This is because, 
with only 3 SO(10) invariant operators, we can redefine the phases of the three 16s 
of fermions to remove the three arbitrary phases. With one more operator, there is 
one additional phase which cannot be removed. A corollary of this observation is 
that this minimal operator set results in just 5 arbitrary parameters in the Yukawa 
matrices of all fermions, 4 magnitudes and one phase. This is the minimal parameter 
set which can be obtained without solving the remaining problems of the fermion 
mass hierarchy, one overall real mixing angle and a CP violating phase. We should 
point out however that the problem of understanding the fermion mass hierarchy and 
mixing has been rephrased as the problem of understanding the hierarchy of scales 
above Mq- Moreover given any particular operator set which fits the low energy 
data we would be obliged at some later time to construct the complete GUT theory, 
including symmetries which forbid additional operators and a consistent description 
of symmetry breaking scales. We leave this problem for future analysis. 

For now we shall describe the detailed analysis of the "22" texture.^ Models with 
"22" texture give the following Yukawa matrices at Mq - 



/ z a C 

z' a C y a Ee<+ x a B | (5) 
V x' a B A 



^As of this writing, we have found no models with "23"' texture which fit the LED. 



with the subscript a = {u, d, e}. The constants x a , x' a , y a , z a , z' a are Clebschs which 
can be determined once the 3 operators ( O23, O22, O12) are specified. Recall, we have 
taken 33 = A 16 3 10 16 3 , which is why the Clebsch in the 33 term is independent of 
a. Finally, combining the Yukawa matrices with the Higgs vevs to find the fermion 
mass matrices we have 6 arbitrary parameters given by A, B,C, E,<f> and tan/3 de- 
scribing 14 observables. We thus obtain 8 predictions. We shall use the best known 
parameters, e,/i, r, c, b, |V^|, as input to fix the 6 unknowns. We then predict the 
values of u, d, s, t, tan/3, \ V c b\, \ V u b\ and J. 

We now show how (within the context of our dynamic principles) we can use the 
LED to guide us towards the theory of fermion masses at Mq. We search for all 
operators ( O23, O22, Ou) with dimension < 10 which fit the data. Using a coarse- 
grained analysis, it is easy to show that just 9 models (with one caveat) may fit the 
data. We now describe this analysis. 

5.1 3rd generation fit - - O33 

Using the values of a and sin 2 9 W evaluated at M z , we obtain a G , M G and a s (M z ). 
There is a theoretical uncertainty in this prediction due to unknown threshold correc- 
tions at both Mq and the weak scale. There is also a 10% experimental uncertainty 
in a s (Mz). Using arbitrary threshold corrections we can obtain any experimentally 
allowed value of a s . We thus allow for all values of a s (Mz) = . 12± .01 self-consistently 
by introducing arbitrary threshold corrections. 
The analysis for the third generation follows — 



where v = 24QGeV and the terms in parentheses are renormalization group factors 
which are implicit functions of both A and a s (Mz). The numerator takes into account 
renormalization from Mz to or m T , and the denominator takes into account the 
running from M G to M z . We use two loop renormalization group equations. 

Using m;, = 4.25±.lGeU and m T = 1.78AlGeV (these are running masses m(m)) 
as input for a given value of a s = .118 (for example) we obtain A from the relation 




(8) 



(6) 



(7) 




Plugging this value of A into the expression for m T we then obtain tan/3. We find 
tan (3 ~ 59 for m b = A.3AGeV. Now using the values of A and tan (3 in the expression 
for m t we find m t (pole) ~ 188GeU. In general, we find values of tan/3 = 54 ± 5 



and m t = 185 ± 15GeV . Both m t and tan/3 increase for increasing values of a s 
or decreasing values of mj. Thus m t = 170GeV and tan/3 = 50 is obtained for 
<x,(Mz) = .110, m b = A.35GeV. Similar results have been obtained previously by 
ALS 6 . Note they find lower values of m t since they allow for values of a s which are 
lower than those presently admissable by the data. 



5.2 2nd generation — O22 

Let us now consider the 2nd generation. We have 4 relations which must be satisfied 
by the LED. 





\v cb \ 




B 

\%u -Ed] ~^ 


~ 1/20 


(9) 


m T 


E 

\y eJ 




X e X e ^ | 


~ 1/17 


(10) 


m s 
m b 


E 
~A 




X d X d ^ | 


~ 1/25 


(11) 


m c 

r^j 


E 
~A 


e i4> 


< B \ 

x u X u ^ | 


~ 10~ 2 


(12) 



We have written these equations using the parameters at Mq, ignoring for the 
moment small renormalization group corrections. Using the first relation we see that 
the ratio B/A ~ 1/10, assuming Clebschs of order 1. The 2nd and 3rd relations 
thus require E/A ~ 1/10. The last relation then requires that the Clebsch, y u « 1. 
Finally, the relation m s = at Mq, first suggested by Georgi-Jarlskog 7 , must be 

incorporated, since it is in good agreement with the LED. We thus conclude that the 
Clebschs, y a , (including RG corrections) should approximately be in the ratio" 

Vu'-Vd- Ve =< 1/3:1: 3. 

The Clebschs, y a , are derived from the operator 022- We have searched over all 
dimension 5 and 6 operators to find solutions to the above Clebsch ratios. We find 6 
solutions — 

16 2 (45 x ) 10 (%^) 16 2 (13) 

16 2 (-j-) 10 (45 B -l) 16 2 
45 x 

16 2 (45 x ) 10 (45 B _ L ) 16 2 

AC 

16 2 10 (^±±) i 6 
V 45 x 

II If the Clebschs are not of order 1 as we assumed, then it is possible that the 4 relations may be 
satisfied with some fine-tuning and a completely different ratio of Clebschs. We have not pursued 
this possibility further. This is our one caveat. 



16 2 10 (45 x 45 B -l) 16 2 
I62 10 16a 

However, note that a// solutions give the same ratio of Clebschs 

Vu ■ Vd ■ Ve = : 1 : 3. 



(14) 



5.3 is£ generation - - Ou 

We can now show that the operator Ou is unique. The first two generations satisfy 
the relations — 



m. 



~ J — tan" /3 — 

m d z d z' d m c \S U J 



\Vcd\ = 




fmd 




l mu e ~i4> 










1 m c 



(15) 
(16) 
(17) 



The first relation is satisfied if 



Z d z 'd ~ Z e z 'e- 



This relation is satisfied if the Clebsch is derived from an SU(5) invariant vev, i.e. 
45x- Since tan/3 is large, the second relation requires 



Z U Z u 



: 3> 



6 or 7 



z d z d 

Finally, the last relation requires 

Zd w z' d . 

The unique operator which satisfies the above 3 relations is 

,45; 

Mp 

5.4 23 



O 



12 



I61 



(18) 



We have now determined, using simple arguments, all but one of the operators. The 
charged fermion mass matrices are given by — 



U 



/ C ON 

C x u B 
V < B A J 



(19) 





( 





-27 C 







D = 




-27 C 


E e** 


x d B 


) 




K 





x' d B 


A 






( 





-27 C 





\ 


E = 




-27 C 


3 E e** 


X e B 






\ 





x' e B 


A 


/ 



(20) 
(21) 

With this form for the mass matrices, the KM element V c b satisfies the relation - 

(22) 



1^1 - X J? 



Si, 



.058 x- 



where x 



m t V Vc S t S 

\ Xu ~ Xd l_ and the last term results from using central values of the input 



parameters. 

Experimentally, we have an upper bound on \V c b\ < .054. We thus require that 
the function of Clebschs, x-, satisfy x < 1- 

We have searched over all dimension 5 and 6 operators for % < 1. We find 9 
possible models with only 3 different values of x 
below. 

X = 2/3 



|, |, or |. The 9 models are listed 



1 

2 
3 
4 



5 
6 



16 2 (45 y ) 10 (-— ) 16 3 
45x 

16 2 (45y) 10 (^f^) 16 s 
45x 

^ (^) 10 (^7) 16 3 



16, 



'45 x 

,45y 



10 16 3 



'45*' v 45 x 
X = 5/6 



4^ 

16 2 (45 y ) 10 (^) 16 3 



'45 



16 2 10 (^) 16 3 



(23) 
(24) 
(25) 
(26) 



'45; 



'45; 



(27) 
(28) 



X = 8/9 



1 



7 16 2 10 (-5-) 16 3 (29) 



8 16 2 10 (-g±) 16 3 (30) 
45 x 



9 16 2 10 (^f^) 2 16 3 (31) 

This is as far as we can get with our coarse-grained search. We must now take 
the 9 distinct models and test the predictions. We are presently in the midst of a 
complete renormalization group analysis, obtaining predictions as a function of the 
input parameters, e, /i, r, c, b, \V c d\ and a s (M z ). For now we present some preliminary 
results. 



6 Preliminary results 

The results in Table 3 are for inputs m c = 1.23GeV, mb = 4.34GeV, |V^| = .221 and 
a s (Mz) = .118 for 4 different models. The model numbers are defined in Eqs. (23 - 
31) 



Table 3. Results for models 3, 6, 8, and 9. 





3 


6 


8 


9 


m d [MeV] 


6.8 


6.9 


6.9 


7.3 


m u /m d 


.82 


.79 


.80 


.79 


m s /m d 


24.8 


24.1 


24.9 


22.4 


m t 


188 


188 


188 


188 


tan (3 


59 


59 


59 


59 


\v cb \ 


.039 


.049 


.052 


.052 




.065 


.064 


.066 


.068 


J x 10 5 


1.8 


2.9 


3.3 


3.7 



Note, | V^b | and J are sensitive to the value of x and so can distinguish between 
the 3 types of models. With better data, \V c b\ will distinguish between these choices. 
At present we can use 6k, the CP violating parameter in K decay, to distinguish 
these models through its dependence on J. We find that models with x — 2/3 tend 
to give too little CP violation. Finally, we see that m s /md and m u /md will also be 
useful in constraining the models. Complete results will be presented in an upcoming 
publication. 1 



7 Conclusion 



We have presented a straightforward method for trying to understand the origin of 
fermion masses. There are four main ingredients — 

• The assumption of a SUSY GUT with its SUSY DESERT implies that the Low 
Energy Data becomes a window into the physics at the GUT scale. 

• The assumption of symmetries (SO(10) plus family symmetries) with a general 
operator analysis allows us to reduce the number of fundamental parameters in 
the theory. 

• With the above assumptions, we obtain predictive theories of fermion masses. 
Given 6 inputs, e,/i,T,c,b, \V c d\, we find 8 predictions for u, d, s, t, tan/3, | V c b\, 
1 14b I an d J- 

• This method allows us to use the LED to systematically find the theory of 
fermion masses at Mq. 

A final note: there is no apriori reason to believe that this scheme should work, 
but if it does work, i.e. if our results agree with the LED, then perhaps we will have 
learned something about the physics at the GUT scale. For example, consider the 
set of possible models 1 - 9 in Eqs. (23 - 31). Each one has the vev 45 x in the 
denominator. This can only occur if this vev gives mass >> Mq to some fermion. 
When this fermion is integrated out in the effective theory at Mq, the vev then appears 
in the denominator of these higher dimension operators. Thus fermion masses and 
mixing angles at low energies may tell us about the hierarchy of symmetry breaking 
scales above M G . 
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